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Abstract 

We introduce the study of neutrino scattering off protons in the deeply virtual kinematics, 
which describes under a unified formalism elastic and deep inelastic neutrino scattering. A 
real final state photon and a recoiling nucleon are detected in the few GeV {\t\ ~ 0.2 — 5 
GeV) region of momentum transfer. This is performed via an extension of the notion of 
deeply virtual Compton scattering, or DVCS, to the case of a neutral current exchange. The 
relevance of this process and of other similar exclusive processes for the study of neutrino 
interactions in neutrino factories for GeV neutrinos is pointed out. 



1 Introduction 



There is no doubt that the study of neutrino masses and of flavour mixing in the leptonic sector will 
play a crucial role for uncovering new physics beyond the Standard Model and to test Unification. 
In fact, the recent discovery of neutrino oscillations in atmospheric and solar neutrinos (see |lj for 
an overview) has raised the puzzle of the origin of the mass hierarchy among the various neutrino 
flavours, a mystery which, at the moment, remains unsolved. The study of the mixing among the 
leptons also raises the possibility of detecting possible sources of CP violation in this sector as well. 
It seems then obvious that the study of these aspects of flavour physics requires the exploration 
of a new energy range for neutrino production and detection beyond the one which is accessible 
at this time. 

For this purpose, several proposals have been presented recently for neutrino factories, where 
a beam, primarily made of muon neutrinos produced at an accelerator facility, is directed to a 
large volume located several hundreds kilometers away at a second facility. The goal of these 
experimental efforts is to uncover various possible patterns of mixings among flavours - using the 
large distance between the points at which neutrinos are produced and detected - in order to study 
in a more detailed and "artiflcial" way the phenomenon of oscillations. Detecting neutrinos at this 
higher energies is an aspect that deserves special attention since several of these experimental 
proposals |ll[2lini require a nominal energy of the neutrino beam in the few GeV region. We recall 
that the incoming neutrino beam, scattering off deuteron or other heavier targets at the detector 
facility, has an energy which covers, in the various proposals, both the resonant, the quasi-elastic 
(in the GeV range) and the deep inelastic region (DIS) at higher energy. In the past, neutrino 
scattering on nucleons has been observed over a wide interval of energy, ranging from few MeV 
up to 100 GeV, and these studies have been of signiflcant help for uncovering the structure of 
the fundamental interactions in the Standard Model. Generally, one envisions contributions to 
the scattering cross section either in the low energy region, such as in neutrino-nucleon elastic 
scattering, or in the deep inelastic scattering (DIS) region. Recent developments in perturbative 
QCD have emphasized that exclusive (see [3j) and inclusive processes can be unifled under a 
general treatment using a factorization approach in a generalized kinematical domain. The study 
of this domain, termed deeply virtual Compton scattering, or DVCS, is an area of investigation of 
wide theoretical interest, with experiments planned in the next few years at JLAB and at DESY. 
The key constructs of the DVCS domain are the non-forward parton distributions, where the term 
non-forward is there to indicate the asymmetry between the initial and final state typical of a true 
Compton process, in this case appearing not through unitarity, such as in DIS, but at amplitude 
level. 

In this work, after a brief summary of the generalities of process, we discuss its generalization 
to the case of neutral currents. We recall that DVCS has been extensively studied in the last few 
years for electromagnetic interactions. The extension of DVCS to the case of neutral currents is 
presented here, while the charged current version will be presented in a forthcoming paper. Also, 
in this work we will just focus on the lowest order contributions to the process, named by us 
Deeply Virtual Neutrino Scattering (DVNS) in order to distinguish it from standard DVCS, while 
a renormalization group analysis of the factorized amplitude, which requires an inclusion of the 
modifications induced by the evolution will also be presented elsewhere. The application of the 
formalism that we develop here also needs a separate study of the isoscalar cross sections together 
with a detailed analysis of the various experimental constraints at neutrino factories in order to 
be applicable at forthcoming experiments. 
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Figure 1: Leading hand-bag diagrams for the process 



2 The Generalized Bjorken Region and DVCS 

Compton scattering has been investigated in the near past by several groups, since the original 
work by Ji and Radyushkin |H1IZ1E|- Previous work on the generalized Bjorken region, which 
includes DVCS and predates the "DVCS period" can be found in [T7| . 

A pictorial description of the process we are going to illustrate is given in Fig[T]where a neutrino 
of momentum I scatters off a nucleon of momentum Pi by an interaction with a neutral current; 
from the final state a photon and a nucleon emerge, of momenta q2 and P2 respectively, while 
the momenta of the final lepton is /'. The process is described in terms of new constructs of the 
parton model termed generalized parton distributions (GPD) or also non-forward (off-forward) 
parton distributions. We recall that the regime for the study of GPD's is characterized by a 
deep virtuality of the exchanged photon in the initial interaction (e + p ^ e + p + ■j) (Q^~2 
GeV^), with the final state photon kept on-shell; large energy of the hadronic system (W'^ > 6 
GeV^) above the resonance domain and small momentum transfers |t| < 1 GeV^. The process 
suffers of a severe Bethe-Heitler (BH) background, with photon emission taking place from the 
lepton. Therefore, in the relevant region, characterized by large and small t, the dominant 
Bethe-Heitler background (~ 1/t) and the 1/Q behaviour of the DVCS scattering amplitude 
render the analysis quite complex. From the experimental viewpoint a dedicated study of the 
interference BH-VCS is required in order to explore the generalized Bjorken region, and this is 
done by measuring asymmetries. Opting for a symmetric choice for the defining momenta, we use 
as independent variables the average of the hadron and gauge bosons momenta 

Pi,2 = P±j qi,2 = qTj (1) 
with —A = P2 — Pi being the momentum transfer. Clearly 

p . A = 0, t = p2 = - ^ (2) 

and M is the nucleon mass. There are two scaling variables which are identified in the process, 
since 3 scalar products can grow large in the generalized Bjorken limit: g^, A ■ g, P ■ q. 

The momentum transfer t = A^ is a small parameter in the process. Momentum asymmetries 
between the initial and the final state nucleon are measured by two scaling parameters, ^ and rj, 
related to ratios of the former invariants 

f = = ^-'^ (3) 

^ 2P-q ^ 2P-q ^ ^ 
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where ^ is a variable of Bjorken type, expressed in terms of average momenta rather than 
nucleon and Z-boson momenta. The standard Bjorken variable x = —qf/{2Pi ■ qi) is trivially 
related to ^ in the t = limit. In the DIS limit (Pi = P2) r] = and x = ^, while in the DVCS 
limit rj = and a; = 2^/(1 + as one can easily deduce from the relations 

We introduce also the inelasticity parameter y = Pi - 1 /{Pi ■ qi) which measures the fraction of 
the total energy that is transferred to the final state photon. Notice also that ^ = ~f measures 
the ratio between the plus component of the momentum transfer and the average momentum. A 
second scaling variable, related to ^ is C = ^^/-Pi^, which coincides with Bjorken x {x = () when 
t = 0. 

^, therefore, parametrizes the large component of the momentum transfer A, which can be 
generically described as 

A = 2^P + A (5) 
where all the components of A are 0{\/t) fT^. 

3 DIS versus DVNS 

In the study of ordinary DIS scattering of neutrinos on nucleons (see Fig. [2j, the relevant current 
correlator is obtained from the T-product of two neutral currents acting on a forward nucleon 
state of momentum Pi 

= uil'h''{-l + Asm^ew + l5)uil) (6) 

where 9w is the Weinberg angle, / and /' are the initial and final-state lepton. The relevant 
correlator is given by 

T,AqI^) = ^ J d'ze^'-^{Pi\T{r,{i)Jm\Pi) (7) 

with u = E — E' being the energy transfered to the nucleon. Pi is the initial-state nucleon 4- 
momenta and gi = / — /' is the momentum transfered. The hadronic tensor W^^ is related to the 
imaginary part of this correlator by the optical theorem. We recall that for an inclusive electroweak 
process mediated by neutral currents the hadronic tensor (for unpolarized scatterings) is identified 
in terms of 3 independent structure functions at leading twist 

Wt^u - \-g^.u + jWi[U,Q ) + Jp l^t^uXaQlPl — — (»j 

where transversality of the current is obvious since Pi^ = Pf — gfPi • qi/q\. 

The analysis at higher twists is far more involved and the total number of structure functions 
appearing is 14 if we include polarization effects. These are fixed by the requirements of Lorenz 
covariance and time reversal invariance, neglecting small CP-violating effects from the CKM ma- 
trix. Their number can be reduced to 8 if current conservation is imposed, which is equivalent to 
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Figure 2: Leading diagram for a generic DIS process 

requiring that contributions proportional to non-vanishing current quark masses can be dropped 
(see also [^). We recall that the DIS limit is performed by the identifications 

MWi{Q\u) = F,{x,Q^) 

vW^{Q\v) = F3(x,Q2), (9) 

in terms of the standard structure functions Fi, F2 and F3. There are various ways to express 
the neutrino-nucleon DIS cross section, either in terms of and the energy transfer, in which 
the scattering angle 9 is integrated over, or as a triple cross section in (Q^, z/, 9), or yet in terms of 
the Bjorken variable x, inelasticity y and the scattering angle {x,y,9). This last case is close to 
the kinematical setup of our study. In this case the differential Born cross section in DIS is given 
by 

^-|ri:-.w^)irH^-. (10) 

the index i denotes the different current contributions, ( i = I7P, \'jZ\, \Z\^ for the neutral 
current) and a denotes the fine structure constant. By 9 we indicate the azimuthal angle of the 
final-state lepton, while y = (Pi ■ qi)/{l ■ Pi) is the inelasticity parameter, and = —qf. The 
factors rji^Q"^) denote the ratios of the corresponding propagator terms to the photon propagator 
squared, 

r^\-\\Q') = 1, 



r2 ' 



^ ^ 2V2na + M| 

(11) 

where Gf is the Fermi constant and Mz is the mass of the Z boson while the leptonic tensor 
has the form 

LU = E [«(fc', A')7m(^?H^?a^75)«(A:, A)]*«(fc',A')7.(^7H^7A^75)«(A:, A). (12) 



4 



In the expression above A and A', denote the initial and final-state helicity of the leptons. The 
indices ii and 12 refer to the sum appearing in eq. (fTn|l 



9'v 
4 



7 
9 a 







1 _L 9 Q^n2 n „Z _ 1 (13) 

2 + ^ sm Uw, 9 A — 2' 

In the case of neutrino/nucleon interaction mediated by the neutral current the dominant 
diagram for this process appears in Fig. [21 A similar diagram, with the obvious modifications, 
describes also charged current exchanges. We recall that the unpolarized cross section is expressed 
in terms of Fi and F2, since F3 disappears in this special case, and in particular, after an integration 
over the scattering angle of the final state neutrino one obtains 



^ ^ 2^S^^M)n^) (2.y^F. + 2(1 - . - . 

where S = 2MEy is the nucleon-neutrino center of mass energy. 

We also recall that in this case the cross section in the parton model is given by 



(14) 



(fa GlME^ 



Ml 



Ml, 



xq\x,Q'')+xq\x,Q''){l-yf 



dxdy 2tc yQ"^ 
where q^{x,Q^) and g°(x, Q^) are linear combinations of parton distributions 



(15) 



+ 



f{x,Q') 



+ 



Uy{x,Q'^) +d^(x,Q^) ^ u{x,Q^) +d{x,Q'^) 



u{x,Q^) + d{x,Q^' 



{Rl + R, 



u.u{x,Q^) + Q^) ^ u{x,Q'^) +d{x,Q^) 



u{x,Q^) + d{x,Q''' 



(16) 



with 



4 2 
1 - - sin^ Ow , Ld= -1 + - sin^ 6w 



4 2 
- sin^ Ow , "^'^ ^ 3 ^^^^ 

(17) 



and we have identified the sea contributions Ug and dg with u and d rispectively. 

Let's now move to the nonforward case. Here, when a real photon is present in the final state, 
the relevant correlator is given by 

T,Ml^)=^ J d'ze^^-{P-^\T{r,{-z/2)j;{z/2)\P+^). (18) 
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Figure 3: The cross section of a neutrino process of DIS-type at x ~ 0.1 with neutral current at 
ultrahigh energy 



The dominant diagrams for this process appears in FigUl We impose Ward identities on both 
indices, which is equivalent to requiring that terms proportional to the quark masses in d ■ Jz are 
neglected. This approximation is analogous to the one performed in the forward case in order to 
reduce the structure functions from 8 to 3 (in the absence of any polarization), imposing symmetric 
trasversality conditions on the weak currents 



p^- — )r^, = o (p'^+ )t,, = o. (19) 



t)^- = » (p^ + - 

The leading twist contribution to DVNS is obtained by performing a collinear expansion of the 
loop momentum of the hand-bag diagram and neglecting terms of order O^A^/Q"^) and M'^/Q'^. 
Transversality is satisfied at this order. Violation of transversality condition in the hand-bag 
approximation is analogous to the DVCS case, where it has been pointed out that one has to 
include systematically "kinematical" twist-3 operators, which appear as total derivatives of twist-2 
operators ^14j in order to restore it. 

For the parametrization of the hand-bag diagram (Fig.[l]) we use the light-cone decomposition 
in terms of 2 four- vectors (n,n), where 



h>' = A(l, 0,0,1) 



= ^(1,0,0,-1) 

= 77,^ = , h ■ n = 1. 



At the same time we set 



= (i+eK + a-Oyn'^-^ 

T/f A^ 

A" = (i + eK + (i + OY^^' + ^ 



q1 = -2^n'^ + ^n^ 

ki" = {k-n) h^ + {k-h) n^' + k'i 

—2 o 

M = (20) 

4 

with = M . We will also use the notation —qf = for the invariant mass of the virtual 
Z boson and we will denote by q the average gauge bosons momenta respectively. 
After a coUinear expansion of the loop momentum we obtain 

74 



(27r 

where we have used the following notations 



n'^^ij T^tr {gr {U. - l') pErgumk)+ 



~gr{Dv-l')rEYgdM'{k)} (21) 



-t. -i 3 

Qu-^e, ga-^e, g ^^^^^^^y 

8 4 

C4 = 1 - - sin^ 9w, D^^l - - sin^ 9w, 



(22) 



{k — aA + gi)^ + ie ' 

^ -g/ +4( (l-g) 
{k-q, + A{l-a)f + ie ^ ' 

where the constant a (a is a free parameter) ranges between and 1. The M matrix is the 
quark density matrix and is given by 

= /Ae*MPlV^?(-«Z/)V'?n(l-«)?/)l^)- (24) 
The index i = u,d runs on flavours. Using a Sudakov decomposition of the internal loop we 



can rewrite 1^ and as 



^ J (27r)4 ^ ^ (27r) 

gdYpDYg {D,-^')M\k)} 

(25) 
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~gY{D,-l')rEY9dM'{k)} (26) 

to which we will refer as the direct and the exchange diagram respectively. It is also convenient 
to introduce two new linear combinations T^'^ = T]^ + = T'^ + T'^ which will turn useful in 
order to separate Vector (V) and axial vector parts (A) of the expansion 

-ggMrib^'por + 7"M1 M"(fc)} + 

(27) 

~ggutv{[YrDl' + 7T£;7l I'M^ik)] + 

~gg,tv {[^^ f^ol' + 7Ti^7l I'K^k)} . (28) 

with Ta including the vector parts {V xV + Ax A) and the axial-vector parts 
{y X A + A X \/). After some algebraic manipulations we finally obtain 



2 / ' " \ 2 



(29) 



where Ci = U^, and 



a(;2) = + , I3[z) = : (30) 

are the "first order" propagators appearing in the factorization of the amplitude. We recall, if not 
obvious, that differently from DIS, DVCS undergoes factorization directly at amplitude level [Til] . 

The parameterizations of the non-forward light cone correlators in terms of GPD's is of the 
form given by Ji at leading twist IH] 
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H{z, e, A')U{P'h^U{P) + E{z, e, A')U{P')^^U{P) + 

H{z,^,A')UiP')r^'UiP) + Eiz,^,A')UiP')^UiP) + (31) 

which have been expanded in terms of functions H, E, H, E j2] and the ellipses are meant 
to denote the higher-twist contributions. It is interesting to observe that the amplitude is still 
described by the same light-cone correlators as in the electromagnetic case (vector, axial vector) 
but now parity is not conserved. 



4 Operatorial analysis 

The operatorial structure of the T-order product of one electroweak current and one electromag- 
netic current is relevant in order to identify the independent amplitudes appearing in the correlator 
at leading twist and the study is presented here. We will identify four operatorial structures. For 
this purpose let's start from the Fourier transform of the correlator of the two currents 

T,, = tJ d'xe''^^P2\T (jJ(a;/2) Jf (-a;/2)) |Pi) , (32) 

where for the neutral and electromagnetic currents we have the following expressions 

/ cos f7vF 

r'\x/2) = M^mY (-^e) Mx/2)+^^{x/2)Y (^e) Mx/2) . 

By simple calculations one obtains 

(P2|T(jJ(x/2)jf (-x/2))|Pi) = 

{PS_uixmguluS{xh,{g!y + g^AY)M-x/2) - 
^a{x/2)gdi,S{x)-f^{g^y + gfja'')^-^/^) + 
fu{x/2)l^.{guv + guAl'')S{-x)gulu^u{x/2) - 
^,{x/2)^,{g^,y + g^,^Y)S{-x)g,^,i;a{xl2)\P,) . 

(34) 

The coefficients used in eqs. fHHl IH^ gy and (7^, are 

guV = \ + \ sin^ dw guA = -\ 

12 1 

gdv = ~2 3 ^^^^ gdA = 2^ (35) 
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and 



9u 2 ) 



9d-^ 



(36) 



are the absolute values of the charges of the up and down quarks in units of the electron charge. 
The function S{x) denotes the free quark propagator 



S{x) 



i:jt 



After some standard identities for the 7's products 
we rewrite the correlators as 



(37) 



(38) 



27r2(x2 -ie)2 



—gd9dV (SnaufiO^ — ie^j,avp0^d^ + QdQdA {S ^avfiO^d^ — ie^au/iO^) 



The x-dependence of the operators in the former equations was suppressed. 
Whence the relevant operators are denoted by 

Of (a;/2, -x/2) = ^(a;/2)7>„(-x/2) + i^,{-x/2)j^M^/2), 
Of (a;/2, -x/2) =ij,{x/2)^'^PM-^m- ^a(-a;/2)7V^a(a;/2), 
Of (a;/2, -x/2) = ^_^(a;/2)7>„(-a;/2) - ^,(-a;/2)7/^^,(a;/2), 
Of (a:/2, -x/2) = ^„(a;/2)7^7^^„(-a;/2) + ^„(-a;/2)7V^a(a;/2) , 

(40) 

where a is a flavour index. 



5 Phases of the Amphtude 

The numerical computation of the cross section requires a prescription for a correct handling of the 
singularities in the integration region at ^ = The best way to proceed is to work out explicitly 
the structure of the factorization formula of the amplitude using the Feynman prescription for 
going around the singularities, thereby isolating a principal value integral (P.V., which is real) 
and an imaginary contribution coming from the d function term. A P.V. integral is expressed in 
terms of "plus" distributions and of logarithmic terms, as illustrated below. The expression of the 
factorization formula of the process in the parton model, in which a{z) and P{z) appear as factors 
in the coefficient functions, is then given by 

Mf,^J^{q,)D{q,)e''*{q,-A) 
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a{z) 



+ 



(5{z)ie^'''''^no.np [^"(z, ^ A')C/(P2)^ I'UiP,) + ^"(z, A^)C/(P2)7' (A • n) U{P^) 
-jOOdDv dz {u^ d}- 



+ 



ggu j^^ dz {-fi^'ny - n^n" + g"") 



a{z) 



H^{z, C, A')U{P2)i l'U{P,) + E^{z, C, A')U{P2) 



i'j^A ■ n 
2M 



U{Pi 



+ 



P{z)i€'""'^han0 



H-{z, e, A')U{P2)i U{Pi) + E^{z, ^, A^)U{P2) 



jSigd J_^dz{u ^ ci}| . 



2M 



(41) 



To handle the singularity on the path of integration in the factorization formula, as we have 
already mentioned, we use the Feynman (ie) prescription, thereby generating imaginary parts. 
In particular, any standard integral containing imaginary parts is then separated into real and 
imaginary contributions as 



/ dz ^^f] = PV C dz^ T ^7rT(±0 

J ZT^±l€ J~l Z^^ ^ ^ 



(42) 



^C^ie J -I z^i 

for a real coefficient T{z). We then rewrite the P.V. integral in terms of "plus" distributions 



where 



1 H{z 
dz — ^ 
-1 z-i 

dzQ{z)H{z) + £ dzQ{z)H{z) + H{i) log 



(43) 



Q{z) = e{-l<z<i)-^ 

= e{-i<z<i) 
Q{z) = e{i<z<i] 



( 0{z<i) 
\{z-i) 



e{-i<z<i) 



z-i)^ 

' {z>0 



{z-0 



(z-O, 



r/ .^ dz 



i). 



(44) 
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and the integrals are discretized using finite elements methods, in order to have high numer- 
ical accuracy. This last point is illustrated in Appendix C, where the computations are done 
analytically on a grid and then the grid spacing is sent to zero. 

We can now proceed and compute the cross section. We define the scalar amplitude 

Mf. = JL,(gi)/^(gi)T^''^e*'^(gi - A) (45) 

where D{qi) is the Zq propagator in the Feynman gauge and Jiepili) the leptonic current 
and we have introduced the polarization vector for the final state photon e'^ . 
In particular, for the squared amplitude we have 

|A^;,,|2 = -L^'^D(gi)'r^.rr (46) 

which is given, more specifically, by 

\Mf = l^^dz l^^dz' {Ki{z,z')a{z)a*{z') + K2{z,z')P{z)P*{z')) (47) 

with Ki and K2 real functions, combinations of the generalized distributions {H, H, E, E) with 
appropriate kinematical factors. Mixed contributions proportional to a{z)[3*{z') and I3{z)a*{z') 
cancel both in their real and imaginary parts and as such do not contribute to the phases. A 
similar result holds also for the pure electromagnetic case. 

After some further manipulations, we finally rewrite the squared amplitude in terms of a P.V. 
contribution plus some additional terms coming from the imaginary parts 

\Mf = P.V. dz dz' {Ki{z, z')a{z)a*iz') + K2{z, z')P{z)P*iz')) 

W (Ki(e, - -0 - + -0) 

W (i^2(e, + ^2(e, -0 + i^2((-e, + i^2(-e, -0) (48) 

which will be analized numerically in the sections below. In order to proceed with the numerical 
result, it is necessary to review the standard construction of the nonforward parton distribution 
functions in terms of the forward distributions, which is the topic of the next section. 

6 Construction of the Input Distributions 

The computation of the cross section proceeds rather straightforwardly, though the construction 
of the initial conditions is more involved compared to the forward (DIS) case. This construction 
has been worked out in several papers [THl [23 [TTl [2HI UHl Uni in the case of standard DVCS, 
using a diagonal input appropriately extended to the non-diagonal kinematics. Different types of 
nonforward parton distribution, all widely used in the numerical implementations have been put 
forward, beside Ji's original distributions, which we will be using in order to construct the initial 
conditions for our process. 

For our purposes it will be useful to introduce Golec-Biernat and Martin's (GBM) distributions 
[TT] at an intermediate step, which are linearly related to Ji's distributions. 
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Figure 4: The relationship between J^'^{X,Q and Ji's function H'^{z,^). 



We recall, at this point, that the quark distributions Hq{z,^) have support in ^ G [—1,1], 
describing both quark and antiquark distributions for z > and z < respectively. In terms of 
GBM distributions, two distinct distributions J-'^{X, () and JF''(X, () with < X < 1 are needed 
in order to cover the same information contained in Ji's distributions using only a positive scaling 
variable (X). In the region X G (C, 1] the functions JF'' and JF^ are independent, but if X < ^ 
they are related to each other, as shown in the (by now standard) plot in Fig. |31 

In this new variable (X) the DGLAP region is described byX>^(|2;|>^), and the ERBL 
region by X < ^ (\z\ < C,)- In the ERBL region, JF'^ and J-'^ are not independent. 

The relation between H{z,C,) and JF'^(X, C) can be obtained explicitly [12] as follows: for 
z G [—(,, 1] we have 



zjH 

1 + e 



and for z G [—1, ^] 



l-C/2 



(49) 



(50) 



1 + l-C/2 ■ 

where i is a flavour index. 

In our calculations we use a simplified model for the GPD's where the dependence can be 
factorized as follows [23 UB] 

H\z,^,A\Q') = Fl{A')q\z,^,Q') 
H\z,^,A',Q') = G\{A')Aq\z,^,Q') 

E\z,^,A',Q') = F^{A')r'{z,^,Q') (51) 

where q\z) and Aq\z) are obtained from the standard non-polarized and longitudinally po- 
larized (forward) quark distributions using a specific diagonal ansatz [3T] . 

The ansatz r*(z, ^) = q\z,^) is also necessary in order for the quark sum rule to hold [H]. 
Analogously, in the case of the distributions [ini 1221 12I| one can use the special model 



m'i — 



X 
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(52) 



valid at small A^, where g^^^ = 1.267, M is the nucleoli mass and is the pion mass, with the 
normalization 



Fm = G1(0) = 1. (53) 

Notice that, analogously to the H distributions, q\z,C,,Q^) and Aq''{z,C,,Q^), which describe 
the = limit of the H-distributions, have support in [—1, 1] and, again, they describe quark 
distributions (for 2; > 0) and antiquark distributions (for z < 0) 

(t{z,^,Q') = -q\-z,^,Q') 

Aq\z,^,Q') = Aq\-z,^,Q'). (54) 

Now we're going to estabilish a connection between the q{z, ^, Q^) and the JF'^(X, Q functions, 
which is done using Radyushkin's nonforward "double distributions" [121 . The construction of the 
input distributions, in correspondence of an input scale Qo, is performed following a standard 
strategy. This consists in generating nonforward double distributions f{x,y) from the forward 
ones (/(x)) using a "profile function" 7r(a;, y) 



f{y,x) =7r{y,x)f{x), (55) 
where we just recall that the vr(?/,x) function can be represented by 

3 [l-\x\r-y' 

"(^'^^ = 4 [1-NP ' ^^^^ 

taken to be of an asymptotic shape (see ref.[Zl[2I|) for quarks and gluons. A more general profile 
is given by 

nix y) - r(25 + 2) [(l-N)^-yT .57. 

^^^'^^ 22''+ir2(6+i) (i-ixD^^'+i 

and normalized so that 



l-\x\ 

dyTc{x,y) = 1. (58) 

-l+\x\ 

b parameterizes the size of the skewing effects starting from the diagonal input. Other choices 
of the profile function are also possible. For instance, the double distributions (DD) defined above 
have to satisfy a symmetry constraint based on hermiticity. This demands that these must be 
symmetric with respect to the exchange y < — > 1 — x — y, and a profile function which respects 
this symmetry constraint is given by ||22j 



6y(l -x-y) 
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Figure 5: GPD's Hu and Hd generated by the diagonal parton distribution with a profile function 
(fBTjl at an initial 0.26 GeV^ 

This symmetry is crucial for establishing proper analytical properties of meson production 
amplitudes. We will be using below this profile and compare the cross section obtained with it 
against the one obtained with ([FTjl . 

Now we are able to generate distributions q{z,^,Q^) in the z variable at = 0, q{z,^,Q^), 
by integrating over the longitudinal fraction of momentum exchange y characteristic of the double 
distributions 



' dx' /' ' dy'5{x' + iy' - z)f{y', x' , Q^). 

-1 J-l+\x'\ 



(60) 



Using (jHUll and the expression of the profile functions introduced above, the GBM distributions 
are generated by the relation 



P^^^X, = y^ dx'n'^ (x', ^(X - x') +x'-l^ q\x') 



(61) 



with a similar expression for the anti-quark distributions in the DGLAP region X > ^ (2; < — ^) 



I r i-c 



T^'\X, = ^J dx'n" Ix', --{X + x') +x' + lj r (|x'| 



(62) 



In the ERBL region, X < ( {\z\ < , after the integration over y, we are left with the sum of 
two integrals 



^''''^(x,c) = - 



rfxV Ix', ^(X - x') + x' - 1 j q%x') 



£ ^ dx'ir'^ (x', ^(X - x') +x' -l\ g^dx'l 
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Figure 6: GPD's flavour singlet combination at 0.26 GeV^ generated with a proflle (f57|) 



^ ^ rfxV (^x', -^{X + x') + x' + 1^ g'^(x') - 
°^ Ix', -^{X + x') + X + 1 j 



(63) 



(64) 



Solving numerically the integrals we obtain the value of the function jF'^s on a grid, and using 
eqs. (|i9|l andEHlwe end up with the numerical form of the H-distributions. We have used diagonal 
parton distribution functions at 0.26 GeV^ [3T| and the results of our numerical implementation 
can be visualized in Figs. El and IHl 

7 The Differential Cross Section 

Our kinematical setup is illustrated in Fig. [HI and we choose momenta in the target frame with 
the following parameterizations 



/ = {E, 0, 0, E) , I' = {E', E' cos 0^ sin 9^, E' sin 0^ sin 9^, E' cos 9^) , 
Pi = (M, 0,0,0) , P2 = {E2, IP2I cos0Arsin6'Ar, IP2I sin0Arsin^7v, I-P2I cos6'iv) 



(65) 



where the incoming neutrino is taken in the positive z-direction and the nucleon is originally at 
rest. The first plane is identified by the momenta of the final state nucleon and of the incoming 
neutrino, while the second plane is spanned by the final state neutrino and the same z-axis. (pu is 
the angle between the x direction and the second plane, while 0Ar is taken between the plane of 
the scattered nucleon and the same x axis. We recall that the general form of a differential cross 
section is given by 



da 



1 



4(/ ■ Pi 



■\Mf^\\27lyS^'\l + Pl-P2-l' 



<fP. 



Q2} 



2ro(27r)3 2P2°(2vr)3 2gO(27r) 



(66) 
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Figure 8: Kinematics of the process v{l)N{Pi) —>■ u{l')N{P2)'-f{q2) 

and it will be useful to express it in terms of standard quantities appearing in a standard DIS 
process such as Bjorken variable x, inelasticity parameter y, the momentum transfer plus some 
additional kinematical variables typical of DVCS such as the asymmetry parameter ^ and A^. 
We will be using the relations 



X 



fl-4 



2x 



2-.(l-A|) 2-. 

(67) 



in the final computation of the cross section. It is also important to note that has to satisfy a 
kinematical constraint 



AL. = -^-;^ 1 + 0^1). (68) 
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One possible cross section to study, in analogy to the DVCS case (22], is the following 



da 



da 



dxdQ^d\A'^\d(f)r dxdyd\A'^\d(f)r 8nQ 



xy^ / 4M^x^\ \ , ^ ,2 



(69) 



where 0,. is the angle between the lepton and the hadron scattering planes. To proceed we also 
need the relations 



where x is given by 



/ ■ n 



(1 + 

2xy 2^ 2 
2xy 



1 - ^'i' 



1 + 



2e 



n ■ gi 



2x 
x-2' 



n- qi = Q'- 



: (2-x) 
8x 



X, 



(70) 



X 



+ 



(71) 



2 4g ' 2 2 

After some manipulations we obtain a simplified expression for similarly to eq. (jTTjl 



1 rl I- 

dz dz Ai{z,z',x,t,Q'^)a{z)a*{z') + A2{z,z',x,t,Q'^)l3{z)p*{z') 

-1 J-l L 



-2Q' (4M' - t) (a; - 2f (x - 1) x'y + {t - AM'^f [x - if a;V 



2M'\Mz' + Q''f {x - 2f [Q'\x - 2f - {AM'' - t){-l + x)^^] V 



-1 



(72) 



where Ai{z, z', x, t, Q^) and A2(-2, -z', x, t, Q^) are functions of the invariants of the process and 
of the entire set of GPD's. Their explicit form is given in Appendix D. As we have already 
mentioned, the (z, z') integration is be done by using the Feynman prescription to extract the 
phases and then using the distributional identities (|i3|l and For numerical accuracy we have 
discretized the final integrals by finite element methods, as shown in Appendix C. 

We have plotted the differential cross section as a function of Q^, for various values of 
and at fixed x values. We have used both the profile given by (f57j) (Figs. l9l fT2|l and the one 
given m eq. (I59ll (Figs. CHI and EI) and their direct comparison in a specific kinematical region 
(Fig. Uni) . The different profiles generate differences in the cross sections especially for larger 
values. Notice also that the DVNS cross section decreases rather sharply with A^, at the same 
time it increases appreciably with x. The results shown are comparable with other cross sections 
evaluated in the quasi-elastic region (~ 10~^ nb) for charged and neutral current interactions, and 
appear to be sizeable. Coherence effects due to neutral current interactions with heavy nuclei, in 
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Figure 9: DVCS cross section at x = 0.1 and center of mass energy ME = 10 GeV^ using the 
profile dsn). 
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Figure 10: DVCS cross section at x = 0.2 and center of mass energy ME = 10 GeV^ using the 
profile (EH). 
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Figure 11: DVCS cross section at x = 0.3 and center of mass energy ME = 10 GeV^ using the 
profile (EH). 
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Figure 14: DVCS cross section at x = 0.3 and center of mass energy ME 
functions generated by the profile function (f59|l . 
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Figure 15: DVCS cross sections at x = 0.2 and center of mass energy ME = 27 GeV^. a using 
profile dSni) 



particular with the neutron component may substantially increase the size of the cross sections, 
with an enhancement proportional to N^, where N is the number of neutrons though an 
accurate quantification of these effects requires a special study |^ which is underway. It is worth 
to emphasize that in the past this contribution had never been included in the study of neutrino- 
nucleon interactions since very little was known about the intermediate energy kinematics in QCD 
from the point of view of factorization. It seems obvious to us that with the new developments now 
taking place in the study of QCD at intermediate energy, especially in the case of the generalized 
Bjorken region, of which the deeply virtual scattering limit is just a special case, it will be of 
wide interest to quantify with accuracy the role of these new contributions for neutrino factories. 
In general, one expects that electromagnetic effects are suppressed compared to the standard 
(hadronic) deeply inelastic cross section, and this has led in the past to a parameterization of the 
intermediate energy cross section as either dominated by the quasi elastic region and/or by the 
DIS region at higher energies, as we have mentioned in our introduction. However, the exclusive 
cross section has some special positive features, one of them being to provide a clean signal for the 
detection of weakly interacting particles, and we expect that this aspect is going to be of relevance 
at experimental level. 

8 Conclusions 

We have presented an extension of the standard DVCS process to the case of one neutral current 
exchange, describing the scattering of a neutrino off a proton in the parton model. We have 
described the leading twist behaviour of the cross section; we have found that this is comparable 
to other typical neutrino cross sections and discussed its forward or DIS limit. We have presented 
a complete formalism for the study of these processes in the parton model. The process is the 
natural generalization of DIS with neutral currents and relies on the notion of Generalized Parton 
Distributions, new constructs in the parton model which have received considerable attention in 
recent years. The possible applications of these new processes are manifold and we hope to return 
in the near future with a discussion of some of the issues not addressed in this work. 
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Appendix A 

The collinear expansion of the internal loop momentum k allows to identify the light cone operators 
appearing in the process at leading twist. We recall at this point that the analysis of the hand-bag 
contribution is carried out exactly as in the electromagnetic case. 

To perform the collinear expansion and isolate the light-cone correlators of DVNS from the 
hand-bag contribution we use the relation 



/ 



d\ dx 



,iX{x—k-n) 



(73) 



inside the expression of T^^ in order to obtain 



r d^k r dX dx 



,iX{x—k-n) 



tr 



r 



27r 
i 



— + + 



+ (1 - a) A -ii + ie 



M{k) 



and therefore 



-I 



tr 



dk ■ n dk ■ n dk"] f dX dx , 

e 



J 2tt J 



r 



1^ + -t"- 



Keeping the leading terms for the loop momenta 



+ (1 - a) A - ^1 + 



r 



k^'-aAi' + q'^ = n''{k-n + 2aC-2C)+n''[k-n-aCM + 



we thus obtain 



-I 



d{k ■ n) d{k • n) dkj_ r dX dx 



{2ny 



J^^^Xix-k.n) 
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tr 



7 



^2 ^.n 7^ + 7^- 



.^r-r2 ,n2 \ 7 



2[k-n- aiM^ + ^) 2 (/t ■ n + (1 - a)eM^ - ^ 



+ 7^ 



2 (A; ■ n + 2< - 20 2 (fc ■ n - 2(1 - a)^ + 2^ 

(P'|V^(-az)^((l-a)z)|P)}. 
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We expand 



/c ■ 2; = {k ■ n) {n ■ z) + {k ■ n) {n ■ z) — k_}^ ■ z_[_ 

and choose a = 1/2. These expansions are introduced in eq. EUand after some manipulations 
the tensor T now becomes 



d{k-n) d{k-h) dkl r dXdx 1 jk-z 

(27r)4 J 2tt 



tr 



+ 7^ 



7 



2 (fc ■ n - f + ^ 



2(A;.n + (l-a)eM'-f^^ 



r 



2 (A; ■ n + 2< - 20 2 (/c ■ n - 2(1 - a)^ + 2^ 

■ {P'\ii{-az)m-a)z)\P)}. 

We also recall that the expansion of the matrix element M.{k) proceeds also in this case as in 
the electromagnetic case 



d'ye''-y{P'\^f{-ay)4'\{l-a)y)\P) = A,^ + A,^S + - (74) 



where the ellipses refer to terms which are of higher twist or disappear in the trace of the 
diagram. 



Appendix B 

The last ingredients needed in the construction of the input distribution functions are the form 
factors Fl and F2. From experimental measurements we know, by a dipole parametrization in the 
small region, that 



GU^') = (1 + K,r'Gl,{A') = ^-'Gl,{A 



1-^1 , G^(A^) = 0, (75) 



where the electric, G^(A^) = F{(A^) + ^^F2(A^), and magnetic form factors (^^^(A^) = 
F{(A^) + F2(A^) are usually parametrized in terms of a cutoff mass my = 0.84 GeV. 

For non-polarized GPD's the valence u and d quark form factors in the proton can be easily 

(P) 



extracted from f)" = 2 + (^'j) Ff and given exclicitely by 
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2F7(A2) = 2Ff (A^) + F7(A2), Ff{A^) = Ff (A^) + 2F^{A^), for 1 = 1,2. 



(76) 



This exploits the fact that proton and neutron form an iso-spin doublet. 
At the scale = 0.9 GeV one can get 



A^^ 



for the valence quarks. For the form factors F we obtain 



where 



1 - 



CA2 



F^ = 



M2 



+ 



F^ 



pd 



M2(l-SA2)2(-l + ^) {l-BA-^Y ' 



CA2 



(1 - BA^f (l - ^) ' 

£;a2 

"m2(1-5A2)2(-1 + ^) ' {l-BA^f ' 



M2 



+ 



(l-SA^)^(l-^)' 



(77) 



(78) 



A = 0.238 B = 1.417 C = 0.447, 
D = 0.835 E = 0.477 F = 0.120. 



(79) 



Appendix C 



In this section we illustrate an analytical computation of the integrals by discretization, using 
finite elements method. We want to approximate with high numerical accuracy integrals of the 
form 



■1 H{z)dz r( ^ H(z) - H{0 , fi ^ H(z) - ff© 



J-1 z — ^ J-1 



e-1 



e + 1 



(80) 



For this purpose we start by choosing a grid on the interval (—1 = xq, , Xn+i — and define 



Ji 



-1 z-i 



E 



^3+1 , H{x 
ax- 



m) 



(81) 
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Performing a simple linear interpolation we get 



71-1 



Xj + l 



^1 = E 

Jxn 



After the integration we are left with 



_Xj-i^l Xj 



+ H{xj+i) 



Xj+l Xj _ 



1 dx 



dx 



dx 

m 

x-i 



(82) 



n—l 

J, = Y.H{x,) 

3=0 

n—l 

i=o 

n—l 
j=0 



-1 + i ^-ti^^ In 



1 + 



^Xjj^l Xj 

^ i-xj 



Xj+l - i 



In 



\Xj+l Xj 

Xj+l - i 



Xj-C 

Xj+l - ^ 



Xj ^ 



Xj-i 

H{xn) + H{C) . 



(83) 



Now, moving to the integral in the interval (^, 1), we introduce a similar grid of equally spaced 
points (^ = I/O, ,yn+i = 1) and define the integral 



As above, after isolating the singularity we obtain 



(84) 



^2 = E 

j=i-'yj 



H{yj) 

+ H{y,) + H{C) 



1 



y-vj 

Vj+i-Vj} 

yi-y 



yi-C 



dy 



+ H{yj+i 

dy T 

e y-^ 



yj+i -yjjj y-i 



dy 



Again, performing the integrations we obtain 



(85) 



J2 = Y.H{yj) 



+ Y.H{yj+i) 

- E^(Oln 



-1 + { y2±^\ In 

Kyj+i-y^J 

l+ ( ^-y^ ] In 



yj+i - C 



\yj+i - yj 
Vj+i - C 



yj-i 
yj+i - ^ 



yj-i 

Collecting our results, at the end we obtain 



yj-i 

+ H{y,) - H{0 . 



H{z)dz 



J1 + J2 + H{C) In 



e + 1 



(86) 



(87) 



25 



We can use the same strategy for the integrals of "+" type defined as follows 



P.V. 



1 H{z)dz f-i H{z) - H{-0 , ,Mz) - Hi-0 



J-i z + ^ "7-1 



dz 



+ / dz- 



e + 1 



This time we call our final integrals Xi and X2. They are given by the expressions 



n— 1 

j=0 
n-l 

j=0 
n-l 

j=0 



In 

\Xj^l Xj 
Xj+i + ^ 



+ { 



In 



Xj +^ 

Xj+i + ^ 



Xj + ^ 



Xj+i 



with a discretization supported in the (—1 = Xq, ,Xn+i — grid, and 



X, = j:H{y,) 

3=1 



-l+ly^^] In 



3=1 

n 

3=1 



( -i- Vj 



Vj+i + C 



In 



on the (— ^ = 1/0, ,2/n+i = 1) grid. As a final result for the "+" integral we get 



•1 H{z)dz 



P.V. I " ^^^7 = + + H{-i) In 



e + 1 



e-1 



Appendix D 



(88) 



(89) 



(90) 



(91) 



In this section we will present the full expression of the functions and A2 which appear in the 
squared amplitude 



Ai{z., z'x., t, 



+ AH',M\4H,{x-l) + E,x')] 

+ Agdgu[{^KHdM^ + AE'^HuA'f + KEdt + E'dEut)x'' 

+ AH'^M\AHa{x - 1) + E^x^) + AH'^M\AH^{x - 1) + E^x^)] 

+ D,U,gdgu[^E'^Edt + 4E^^„t - AE'^E^tx - AE'^Ejx + AE'^E^M^x^ 

+ AE'^E^M^x"" + AE'^HdM^x^ + AE'^H^M^x^ + K^rf^a:' + ^^^.to;' 

+ AH'^M\AHa{x - 1) + E^^x^) + AH'^M\AH^{x - 1) + £;„x2)] 
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+ gl[4ElEjU^ - AE'^E^tUlx + WKHuM^x^ + AE^Ejx'' 
+ AE'^EuM^Ulx^ + AE'^H^M^Uy + E'^EjUy + IQE'^M^AH^^x - 1) + 
+ AH'^M^Ul{AHu{x - 1) + E^x^)] 

+ Dlgl[AH'^M\m,{x - 1) + E,x^) + E'^iAH^M^x^ + E,(t [x - 2f + 4M2x2))]' 

(92) 

and for Aoi t) we get a similar result 



A2{z, z', X, t, = 4^4g2 ^^^^^[4^^^^^ + ^E'^EJ - 16D,H'^HdM^U, - WD^H'^H^M^U, - AE'^E^tx 

- AE'^Eutx + WDyHuHdM^U^x + IQDyH'^HuM^U^x + AE'^EaM^x^ + AE'^E^M^x'^ 

+ AE'^HdM^x^ + AE'^HuM'^x^ + E'^Edtx^ + E'^Eutx^ + AD^E^H'^M'^U^ 

+ AD^EdKM^U^x'^ + AD^E'^HdM^U^x^ + AD^E'^HuM'^Uy + D^E'^EatUy 

+ D.E'dEutUy + AHiM\AHd{x - 1) + E^x') + AH'dM\AH^{x - 1) + E^x')] 

+ gl[AH'^M\AHd{x - 1) + E^x^) + Dl{E',{AHdM^ + Erft)^^ 

+ AH'dM^{AHd{x - 1) + E^x')) + E'diAHdM^x'' + Ed{t{x - 2f + 4M2x2))] 

+ AH'^M\AH.,{x - 1) + E„x2)) + E'^{AH.,M^x^ + E„(t(x - 2)^ + 4M2x2))]" . 

(93) 
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